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ABSTRACT 

The  rotor  system  of  the  proposed  XH-2/CCR  (Circulation 
Control  Rotor)  prototype  aircraft  and  the  state  variable 
format  of  the  airframe  equations  of  motion  are  described. 
Through  a  study  of  the  eigenvalues  and  eigenvectors  of  the 
basic  airframe,  the  effects  of  uncoupling  and  cross-coupling 
the  helicopter  equations  of  motion  were  analyzed.   The 
control  matrix  for  the  CCR  was  then  generated  using  an 
applicable  pneumatic  lead  angle  for  Coanda  blowing.   Feed- 
back gains,  to  give  the  aircraft  acceptable  flying  qualities, 
were  calculated  for  the  Stability  Augmentation  System  (SAS) 
and  finally  the  aircraft's  sensitivity  to  changes  in  pneu- 
matic lead  angles  were  studied.   The  programs  used  in  the 
above  study  (capable  of  handling  up  to  10  x 10  matrices) 
are  listed  for  both  the  IBM-360  and  HP-9830  mini-computer. 
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I.   INTRODUCTION 

In  the  early  1970* s  the  David  W.  Taylor  Naval  Ship 
Research  and  Development  Center  (DTNSRDC)  initiated  a 
research  program  into  the  feasibility  of  incorporating  a 
Circulation  Controlled  Rotor  (CCR)  system  in  a  Navy  heli- 
copter.  The  concept  of  the  CCR  and  of  improving  airfoil 
lift-to-drag  ratios  using  Coanda  flows  had  been  proven 
earlier  by  some  of  the  world's  leading  aerodynamicist' s 
and  as  early  as  1959  Dorand  had  published  in  the  Journal  of 
the  Helicopter  Association  of  Great  Britain  an  article  on 
the  application  of  a  jet  flap  to  control  a  helicopter  rotor 
(Ref.  1).   Studies  and  tests  continued  throughout  the  1960's 
with  more  papers  published  in  both  the  United  States  and 
in  Europe  on  the  improved  performance  and  possible 
applications  of  a  CCR  (Ref.  2) . 

The  Aviation  and  Surface  Effects  Division  of  DTNSRDC 
continued  the  research  with  further  tests  involving  detailed 
pressure  measurements  of  two-dimensional  elliptical  sections 
in  their  15  x 20-inch  subsonic  wind  tunnel.   These  tests 
reconfirmed  that  extremely  high  lift-to-drag  ratios  could 
be  achieved  by  tangentially  injecting  air  through  a  slot 
in  the  trailing  edge  of  an  airfoil.   The  results  of  these 
tests  were  incorporated  in  both  two- and  four- b laded  model 
rotor  systems  for  evaluation  in  the  DTNSRDC  8x10  foot 
windtunnel  (Ref.  3)  . 


The  incorporation  of  a  CCR  system  in  a  full-size 
aircraft  could  conceivably  offer  other  advantages  over  the 
conventional  rotor  system,  in  addition  to  the  potentially 
improved  aerodynamic  performance  traits.   The  conventional 
rotating  mechanical  swashplate  system  would  be  replaced 
by  a  non-rotating  pneumatic  plate  in  a  plenum  chamber 
located  in  the  blade  hub  region.   Collective  control  would 
be  accomplished  by  changing  the  plenum  chamber  pressure, 
which  increases  or  decreases  the  Coanda  blowing  equally  at 
all  blades  via  the  individual  blade  supply  or  collector 
tubes.   Cyclic  control  would  be  provided  by  tilting  the 
pneumatic  swashplate  so  that  there  is  an  azimuthal  variation 
in  Coanda  blowing  in  each  blade.   This  variation  in  blowing 
is  a  result  of  the  changes  in  volume  of  air  allowed  to  the 
collector  tubes  because  of  changes  in  the  gap  between  the 
collector  tubes  and  the  non-rotating  swashplate.   This 
non-rotating  swashplate  and  variations  in  Coanda  blowing 
would  replace  the  mechanical  cyclic  feathering  required  by 
conventional  rotors  and  therefore  eliminate  the  vibrations 
caused  by  this  one-per-revolution  cyclic  mechanical  move- 
ment.  Another  important  point  to  recognize  is  that  the 
CCR  concept,  with  the  Coanda  blowing,  will  dictate  a  tor- 
sionally  stiff  rotor  blade  or  rigid  rotor  system.   This 
is  a  result  of  the  disparity  between  the  two  lift  generation 
centroids.   The  center  of  pressure  due  to  Coanda  blowing 
is  near  the  blade  midchord  region,  while  the  blade  aerodynamic 
center  remains  near  the  rotor  blade  quarter  chord  point. 
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The  proposed  helicopter,  with  a  simple  hub  and  rigid 
rotor  blades  free  of  flapping  and  lag  hinges,  would  result 
in  a  relatively  clean  aerodynamic  hub  system.   The  reduction 
in  rotor  and  hub  drag  would  be  beneficial  to  the  helicopter 
from  a  performance  standpoint.   The  reduction  in  moving 
parts  in  the  hub  and  blade  system  would  also  mean  a  quieter 
helicopter  with  a  lower  vibration  level  than  that  of  a 
conventional  rotor  system.   This  latter  effect  has  a  favor- 
able potential  of  improving  the  "ilities"  (maintainability 
and  reliability)  for  helicopter  operations. 
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II.   BACKGROUND 

Early  in  1973  the  Naval  Air  Systems  Command  (NAVAIR) 
contracted  the  Lockheed  Aircraft  Company  and  Kaman  Aero- 
space Corporation  to  investigate  the  feasibility  of  devel- 
oping a  full-scale  flightworthy  Circulation  Controlled 
Rotor  demonstrator  aircraft.   In  the  summer  of  1974,  some 
twelve  months  later,  both  companies  returned  reports  to 
DTNSRDC  and  NAVAIR  stating  that:   "the  concept,  while 
innovative,  is  completely  safe  in  operation"  (Ref.  4)  and 
"that  there  is  no  fundamental  flaw  or  deficiency  in  the 
CCR  concept  and  that  construction  of  a  full  scale  CCR 
helicopter  is  feasible  and  practical"  (Ref.  5) .   Lockheed 
Aircraft  proposed  the  use  of  its  L2  86/CCR  while  Kaman 
suggested  "that  the  Kaman/Navy  H-2  aircraft  is  an  ideal  test 
vehicle  for  the  CCR  concept"  (Ref.  4) . 

With  the  additional  goal  of  being  able  to  incorporate 
a  CCR  system  on  any  off-the-shelf  helicopter  with  no  major 
airframe  or  equipment  changes ,  NAVAIR  awarded  a  contract 
to  Kaman  Aerospace  to  "develop,  build  and  test"  a  prototype 
CCR  vehicle  incorporating  the  use  of  the  Navy/Kaman  H-2 
aircraft.   This  technology  demonstration  aircraft  will 
tentatively  be  designated  as  the  Navy  XH-2/CCR. 

Preliminary  studies  by  Kaman  promoted  the  belief  that 
acceptable  flying  qualities  would  be  sustained  with  the 
installed  Stability  Augmentation  System  (SAS)  of  the  Kaman 
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H-2  with  only  minor  changes  in  the  gains  of  the  feedback 
amplifiers  (Ref .  5) .   Acceptable  flying  qualities  does  not 
necessarily  mean  all  stable  roots  of  the  aircraft  motion 
modes,  since  a  weak  oscillatory  instability  with  a  time- 
to-double  amplitude  of  greater  than  three  seconds  can  be 
tolerated  by  a  proficient  rotor-wing  aircraft  pilot.   The 
objectives  of  this  research  was  to  confirm  that  Kaman's 
beliefs  were  in  fact  true  and  to  find  a  suitable  feedback 
law  for  the  SAS  of  the  XH-2/CCR  such  that  within  the  air- 
craft's flight  envelope  the  aircraft  flying  qualities  will 
be  acceptable  to  the  evaluation  pilot. 
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III.   DISCUSSION 

The  study  of  the  helicopter  flight   dynamics  were  con- 
ducted using  the  conventional  non-dimensionalized  state 
variable  format  of  the  aircraft  linearized  equations  of 
motion  (Ref .  6) ,  modified  to  allow  coupling  of  the  longi- 
tudinal with  the  lateral-directional  motions.   This  modi- 
fication is  a  fairly  elementary  record-keeping  operation 
when  using  state  vector  formulations.   The  basic  plant 
matrix,  A,  of  the  aircraft  linearized  equations  of  motion 
in  the  state  variable  format  is  given  in  Appendix  A. 

The  stability  derivatives  for  the  XH-2/CCR  airframe 
were  computer  generated  by  the  contractor  using  the 
MOSTAB-HFA  program  (Ref.  7)  modified  for  the  pertinent 
characteristics  of  the  SH-2F  airframe  and  the  XH-2/CCR 
main  rotor  system.   The  flight  conditions  analyzed  were  for 
l.Og  level  flight  at  sea  level  standard  conditions.   The 
aircraft  gross  weight  was  given  as  11,000  pounds  and  a 
rotor  tip  speed  of  615  feet  per  second  (2  67  RPM)  was  used 
throughout  the  calculations.   Stability  derivatives  were 
generated  for  airspeeds  of:   Zero  (hover),  35,  72,  110, 
and  130  knots.   These  derivatives  were  computed  in  May  of 
1976  and  then  updated  in  November  of  the  same  year.   The 
calculations  made  in  this  research  effort  are  based  on 
the  updated,  November  1976,  data. 
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A.   THE  PLANT  MATRIX 

The  plant  matrix,  A,  was  developed  for  the  longitudinal 
and  lateral-directional  components  and  then  the  fully- 
coupled  equations  of  motion  using  the  contractor  generated 
stability  derivatives.   The  plant  matrix,  A,  was  partitioned 
into: 


A  = 


An   ! 


L21   I 
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22 


where   [A, ^]   represented  the  coefficients  of  the  long- 
itudinal stability  derivatives  and   tA22^   represented 
the  coefficients  of  the  lateral-directional  stability 
derivatives.   The  cross-coupling  stability  derivatives 
were  represented  by  the  coefficients  of   [A-|->]   an<^   ^Aoi^ 

The  homogeneous  form  of  the  state  equations  took  on 
the  conventional  form  of : 


r  . 

X. 


X 


-J— 


.  2J 


=   A  x   =   A^ 


.\ 


where  the  state  vector,  x,  included  the  partitioned 
longitudinal  airframe  state  vector: 


x1   =   [u,  w,  q,  0] 
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and  the  lateral-directional  airframe  state  vector 
respectively: 


T 
x0   =   [p,  r,  v,  <J>] 


Normally  in  airframe  dynamics  the  fact  that  most  "fixed 
wing"  aircraft  approach  symmetric  conditions  enables  one 
to  eliminate  the  cross  coupling  terms   [A.J   and   [A-. ] 
and  then  to  represent  the  aircraft  by  the  decoupled  equations 
of  motion  in  the  longitudinal  and  lateral-directional  modes. 
Neglecting  these  cross-coupling  terms  allows  analysis  of 
the  individual  fourth-order  systems  and  the  size  and  shape 
of  the  applicable  longitudinal  and  lateral-directions  modes 
changes  only  slightly  when  the  more  complicated  computations 
are  made  for  the  fully  coupled  eighth-order  systems.   Unfor- 
tunately, helicopters  do  not  enjoy  these  conditions  of 
symmetry  and  the  effects  of  completely  cross-coupling  the 
longitudinal  and  lateral-directional  equations  are  quite 
significant. 

The  fact  that  the  helicopter  had  large  cross-coupling 
effects  was  shown  by  the  significant  change  in  the  eigen- 
values for  the  uncoupled  and  fully-coupled  computations 
and  modal  identification  could  not  be  completed  from  the 
results  of  the  cross-coupled  forms  alone.   The  conclusion, 
therefore,  was  that  the  stability  problem  could  not  be 
solved  by  studying  only  the  uncoupled  components  of  the 
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equations  of  motion  but  must  be  undertaken  using  the 
fully  coupled,  eighth-order  equations  of  motion. 

Using  the  Basic  Matrix  Control  Theory  (BASMAT)  com- 
puter program  (Ref.  8)  modified  for  the  HP-9830  mini- 
computer and  the  IBM-360  digital  computer  located  at  the 
Naval  Postgraduate  School's  W.  R.  Church  Computer  Center; 
calculations  were  made  at  the  above  stated  airspeeds  and 
the  eigenvalues  and  eigenvectors  of  the  plant  matrix,  A, 
were  obtained  for  both  the  uncoupled  and  fully-coupled 
equations  of  motion.   The  revised  programs  for  both  the 
IBM-360  and  the  HP-98  30  are  included  as  part  of  this  paper 
in  Appendix  B  and  C  respectively. 

The  existence  of  a  longitudinal  unstable  root  in  the 
uncoupled  form  was  confirmed  at  all  calculated  speeds  and 
the  identification  of  the  different  modes  in  the  fully- 
coupled  eighth-order  system  was  completed  by  slowly  intro- 
ducing the  cross-coupling  derivatives  into  the  aircraft 
equations  of  motion.   The  eigenvalues  and  eigenvectors  were 
traced  and  modes  identified  by  letting  A,,,   and  A..  _ 
equal  zero  in  the  eighth-order  system  and  then  slowly 
increasing  their  values  until  they  reached  the  final  values 
of  the  basic  plant  matrix,  A.   The  values  of  the  uncoupled 
eigenvalues  from  the  fourth-order  solutions  and  the  final 
eighth-order  system  results  are  included  in  Table  I. 

A  modified  root-locus  is  shown  in  Figures  1  and  2 
for  the  airspeed  conditions  of  hover  and  130  knots  respec- 
tively.  The  trajectory  of  the  roots  is  shown  as  the  amount 


17 


TABLE    I 


Note: 


Hover 


Eigenvalues   of   Basic   XH-2/CCR  Airframe 

Uncoupled  airframe  listing  shows  four  longitudinal  roots  first 
followed  by  four  lateral-directional  roots. 


(uncoupled) 

i  0.10077 


X  =     0.00159  ± 
=  -0.23713 
=  -2.64654 

=     0.02949  ±  i  0.27626 
=  -0.44778 
=  -7.71637 

35  Knots     (uncoupled) 

A  =  -0.00971  ±  i  0.16470 
=  -0.30077 
=  -2.37101 

=  -0.24771  ±  i  0.81760 
=  -0.05339 
=  -6.61402 

72  Knots     (uncoupled) 


X  = 


-0.79763 

-0.14599 

-0.11719 

-2.19846 

-0.41632  ± 

-0.05296 

-6.10178 


i  1.42319 


X  = 


(coupled) 


X  = 

-0.28206  ± 

i  0.01601 

= 

0.01867 

= 

-0.05189 

= 

-0.00634  ± 

i  0.20903 

-5.18781  ± 
(coupled) 

i  2.53260 

X  = 

0.01046  ± 

i  0.13193 

= 

-0.24564 

= 

-4.50717  ± 

i  2.34599 

= 

-0.27651  ± 

i  0.85642 

= 

-0.06194 

(coupled) 

-0.05862  ± 
-0.27382 
-4.40993  ± 
-0.42956  ± 
-0.05781 


i  0.15912 

i  2.33811 
i  1.53113 


110  Knots   (uncoupled) 

^  =  -0.38109  ±  i  0.26126 
=     0.31503 
=  -2.92894 

=  -0.54931  ±  i  1.97156 
=  -0.04699 
=  -5.59397 


X  = 


(coupled) 

-0.15433 
0.36151 
-4.52643 
-0.53360 
-0.04816 


±  i  0.25617 

±  i  2.93431 
±  i  2.10332 


130  Knots    (uncoupled) 

X  =  -0.30712  ±  i  0.29749 
=     0.31242 
=  -3.02559 

=  -0.63733  ±  i  2.22591 
=  -0.05636 
=  -5.00383 


X  = 


(coupled) 

-0.07151 
0.41211 
-4.35708 
-0.56744 
-0.08227 


±  i  0.20758 

±  i  2.81962 
±  i  2.34661 
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of  cross-coupling  increased  from  0.0  to  1.0  with  the  latter 
limit  corresponding  to  a  fully  cross-coupled  system.   The 
Dutch-roll  and  long  period  oscillatory  roots  are  omitted 
for  sake  of  clarity,  but  their  values  do  not  vary 
significantly  with  cross-coupling  as  shown  in  Table  I. 

The  uncoupled  non-oscillatory  lateral-directional  roots 
may  be  identified  as  spiral  and  roll  subsidence  roots  by 
both  the  mode  shape  and  time  constants  using  familiar 
analogies  from  fixed  wing  aircraft.   The  remaining  two 
non-oscillatory  roots  from  the  uncoupled  longitudinal 
degrees  of  freedom  would  normally  correspond  to  a  short 
period  situation  in  fixed  wing  aircraft. 

It  was  observed  that  the  low  time  constant-real  roots, 
one  each  from  the  longitudinal  and  lateral-directional 
degrees  of  freedom  respectively  remain  almost  invariant 
with  the  amount  of  cross-coupling  until  they  reach  the 
neighborhood  of  full  cross-coupling.   At  that  time  the 
longitudinal  root  becomes  weakly  unstable  with  a  time-to- 
double  amplitude  of  approximately  37  and  1.68  seconds  at 
hover  and  130  knots  respectively.   The  latter  situation 
definitely  required  improvement  by  means  of  stability 
augmentation.   These  two  real  roots  could  have  been  expected 
to  coalesce  into  an  oscillatory  pair  as  cross-coupling 
varied,  but  possibly  the  close  proximity  to  the  almost 
invariant  oscillatory  long  period  and  Dutch-roll  roots 
prevented  this  action  from  occurring. 
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The  second  real-longitudinal  root  and  the  roll  subsi- 
dence root  may  be  observed  to  coalesce  into  a  pair  of 
complex  conjugate  (oscillatory  roots)  almost  mid-range  in 
the  cross-coupling.   The  coalesced  roots  have  a  mode  shape 
similar  to  a  pendulum  type  of  motion,  but  it  was  noted  that 
this  mode  (in  the  fully-coupled  situations)  was  quite 
heavily  damped. 

The  aircraft  mode  shapes  are  defined  in  Table  II  for 
hover  and  130  knots  velocities  in  both  the  uncoupled  and 
fully  coupled  situations.   Without  these  mode  shapes, 
there  would  be  difficulties  involved  in  interpreting  the 
Characteristic  root  migrations  as  shown  in  Figures  1  and 
2.   The  symbol  (O)  is  used  on  the  figures  to  identify  a 
longitudinal  plunging  subsidence  mode  which  will  be  spotted 
in  Table  II  with  the  (w)  velocity  perturbation  being  the 
dominant  term.   This  root,  which  in  conventional  airframe 
systems  would  be  combined  with  the  rocking  mode  root  (  □  ) 
to  yield  the  conventional  oscillatory  short  period  mode, 
is  the  "culprit"  which  becomes  unstable  in  the  fully  coupled 
situation.   As  will  be  noted  in  Figures  1  and  2,  the  longi- 
tudinal rocking  and  the  lateral  roll  subsidence  modes 
.coalesce  into  an  oscillatory  pendulum  type  mode. 

B.   THE  CONTROL  MATRIX 

In  addition  to  the  definition  of  the  airframe  plant, 
the  contractor  provided  the  MOSTAB  generated  control 
matrices.   The  derivative  information,  presumably  the  result 
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TABLE    II 
MODE    SHAPE    SUMMARY 


A.     Hover  Conditions 

1.     Zero  Cross-Coaling 

a.  Longitudinal-oscillatory  long  period 


A=     0.0016    ±i  0.1008 

r 


Period  =62.3  sec 


1.0000;  arg.  0.0  deg. 

0.0196;  arg.  -116.8  deg. 

0.0003;  arg.  83.8  deg. 

0.0031;  arg.  -  5.4  deg. 


b.  Longitudinal-plunging  subsidence 


T2  =  435.94  sec. 


A  =  -0.2371 

0.0700 
1.0000 
0.0000 
0.0000 


T,  =  2.92  sec. 


J 


c.  Longitudinal- fore/aft  rocking  subsidence 


A  =  -2.647 

1.0000 
-0.0090 
-0.1830 

0.0690 


^ 


T..  =  0.26  sec. 


I 


J 


d.  Lateral-Directional-Dutch  Roll 

X  =  0.0295  ±  i  0.2762      Period  =  22.8  sec. 

r 


T2  =  23.4  sec. 


0.0025;  arg.  166.4  deg. 

0.0250;  arg.  -  34.4  deg. I 

1.0000;  arg.  0.0  deg. f 

0.0090;  arg.  82.4  deg. 


J 


e.  Directional-Spiral  subsidence 


A  =  -0.448 


L 


r 0.0044^ 


H 


T-  =  1.55  sec. 


-0.1650 
1.000 
0.0099 


f. 


Lateral-Roll  Subsidence 
X  =  -7.716 


S 


P 
v 


Ut 


0000 
,0294 
0.8669 
-0.1296 


T^  =  0.090  sec. 


j 
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TABLE   II    (Continued) 


A.     2.     Full  Cross-coupling 

a.     longitudinal-oscillatory  long  period 


A  =  -0.2821  ±  i  0.0160 


r  ^ 

u 

w 

q 


< 


p 

r 
v 


v-  j 


1.0000;  arg 

0.7102;  arg 

0.0025;  arg 

,_  J  0.0089;  arg 
^0.0008; 
0.0490; 

0.1924;  arg 

0.0028;  arg 


arg. 
arg. 


0.0  deg. 
14.1  deg. 

173.4  deg. 

-  3.3  deg. 
176.0  deg. 

174.5  deg. 

182.6  deg. 

-  0.7  deg. 


Period  =  392.7  sec. 
^  T,  =  2.46  sec 


J 


b.  Cross-coupled  Pendulum  Mode 
A  =  -5.1878  ±  i  2.5326 


r  "N 
u 

w 


<lt=< 


p 

r 
v 


0.6181;  arg, 

0.0798;  arg. 

0.4219;  arg. 

0.0731;  arg, 

0.7618;  arg, 

0.0386;  arg, 

1.0000;  arg. 

0.1320;  arg. 


-141.88  deg. 

155.40  deg. 

2.14  deg. 

-155.84  deg. 

-  42.92  deg. f 

177.95  deg. 

0.0  deg. 

163.10  deg. 


Period  =2.48  sec. 
-\  T,  =0.13  sec 


j 


c.  Lateral-Direction  -  Coupled  Dutch  Roll  Mode 


A  =  -0.0063  ±  i  0.2090 


Period  =  30.1  sec. 


r  *s 
u 

w 

q 


< 


r. 


1.0000; 
0.0162; 
0.0014; 


arg. 
arg. 
arg. 


P 
r 

v 


_)  0.0065;  arg. 
f   l 0.0006;  arg. 
0.0126;  arg. 
0.4099;  arg, 
0.0028;  arg. 


0.00  deg. 

-171.87  deg. 

3.04  deg. 

-  88.70  deg. 

145.88  deg. 

1.69  deg. 

33.20  deg. 

122.38  deg. 


"^  T 


109.26  sec. 


>- 


d.  Lateral- 
A  =  -0.0519 


< 


u 

w 

q 

9 

p 

r 
v 


H 


j 


1.0000^ 
-0.0358 
-0.0000 

0.0016 

-o.oooo  r 

-0.0077 

-0.1992 

0.0003 
v.      J 


T,  =  13.36  sec. 
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TABLE    II    (Continued) 


e.     longitudinal 

X  =     0.0187 

^1.0000^ 

-0.0224 

-0.0000 

,"0.0006 

«  1-0.0000 
-0.2016 
-0.0001 


T2  =  37.13  sec. 


< 


u 

w 

q 

e 

p 

r 


> 


J 


B. 


130  Knots 
1. 


Zero  Cross-coupling 

a.     Longitudinal-oscillatory  Long  Period  Mode 

X  ='-0.3071  ±  i  0.2975       Period  =  21.1  sec. 


T,  =  2.26  sec. 


0.5812;  arg. 

1.0000;  arg. 

0.0026;  arg. 

0.0061;  arg. 


b.  Longitudinal  - 


-67.70  deg. 

0.0  deg. 

39.45  deg. 

-96.46  deg. 


> 


j 


X  = 


0.3124 

l.oocxf) 

-0.7721  [ 
-0.0043 
0.0138 


T~  =  2.22  sec. 


f 


c. 


J 

Longitudinal  - 

X  =  -3.0256 

0.0185*") 
1.0000  I 
-0.0108  f 
0.0036 

Lateral-Directional  Dutch  Boll  Mode 


T1  =   0.23  sec. 
2 


X  =  -0.6373  ± 


^ 


< 


r 
v 


>=< 


0.0185;  arg. 
1.0000;  arg. 
0.0107; 
0.0036; 


i  2.2259   Period  =2.8  sec. 

0.02  deg?) 

0.0  deg. I 
arg.  179.96  deg.  ' 
arg.  -     0.0     deg. 


T1  =  1.09  sec. 
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TABLE    II    (Continued) 


e.     Lateral- 


-0.0077) 

J   0.0203  1 

1.0000  r 

0.1363 


V 


f.     Lateral- 
X  =  -5.0038 


r 


1.0000 
-0.0381 
-0.0674  ( 
-0.1999 


> 


j 


B.     2.     Full  Cross-coupling 


T  =  12.29  sec. 


T,  =  0.14  sec. 


a.  Longitudinal-oscillatory  Long  Period  Mode 

X  =  -0.0715  ±  i  0.2076   Period  =  30.3  sec.  T,  =  9.69  sec 


u 

w 


< 


arg 
arg 


P 
r 
v 

.  <J>  - 


>-< 


1.0000; 

0.4987; 

0.0014;  arg 

0.0064;  arg 

0.0045; 

0.0027; 

0.3260; 


\_ 


0.2051; 


arg 
arg 
arg 
arg 


0.0  deg. 

10.05  deg. 

34.18  deg. 
-  74.83  deg. 
167.65  deg. 

81.22  deg. 

25.33  deg. 
883.34  deg. 


-\ 


> 


b.  Lateral-Directional  -  Coupled  Dutch  Poll  Mode 

X  =   -0.5674  ±  i  2.3466    Period  =  2.7  sec.  T,  =  1.22  sec 


r  ^ 

u 
w 

q 


r0.0111;  arg.  -109.36  deg.^ 
0.2045;  arg.  -  20.66  deg. 
0.0023;  arg.  -  74.39  deg. 


^H 


P 
r 
v 


0.0009;  arg. 

0.0045;  arg. 

0.0110;  arg. 

1.0000;  arg. 

0.0019;  arg.   56.22  deg. 


88.99  deg. 

159.81  deg. 

-  76.86  deg. 

0.0  deg. 


J 
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TABLE    II    (Continued) 


c.      Cross-couple  Pendulum  Mode 

X  =  -4.3571  ±  i  2.8196         Period  =  2.2  sec. 


u 
w 

q 


< 


p 

r 
v 


M 


L 


0.0177 
1.0000 
0.0214 
0.0041 
0.0412 
0.0066 
0.2695 
0.0079 


arg, 
arg, 
arg. 
arg. 
arg. 
arg. 
arg, 
arg. 


12.86  deg. 

0.0     deg. 

-66.01  deg. 

-  3.95  deg. 

82.02  deg. 
-8.08  deg. 

37.04  deg. 
-65.07  deg. 


^ 


T.   =0.16  sec. 


> 


d.  Longitudinal- 
X  =  0.4121 


T-  =  1.68  sec. 


1T> 

w 

q 


< 


p 

r 
v 

<$> 


>  A. 


-0.7732 
1.0000 
0.0006 
0.0014 
-0.0031 
-0.0014 
-0.0342 
-0.0075 


*> 


> 


J 


J 

d.  Lateral 


X  =  -0;0823  T,  =  8.42  sec. 

AT}  C  1.000(T|2 
w  0.0781 
-0.0000 
0.0009 
-0.0041 
r  0.0076 
v  0.5428 
4>      0.0504 


H 


> 
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of  pneumo-dynamic  modeling  in  the  hub  and  blade  blowing 
sections ,  was  supplied  as  airframe  forces  and  moments  per 
unit  (3060  psf)  pressure  variation  at  the  pneumatic  swash- 
plate  referenced  to  the  aircraft  axis.   The  cyclical  varia- 
tion of  the  plenum  pressure  upon  control  was  given  by: 


C   (2)   .  .  The  coefficient  of  the  cosine  4;type 
variation  in  the  control  matrix. 

C   (3)   ..  The  coefficient  of  the  sine  ip  type 
variation  in  the  control  matrix. 


where,  ip,  represents  the  blade  azimuthal  angle.   One  can 
visualize  the  longitudinal  and  lateral  cyclic  controls  as 
rotating  the  pneumatic  swashplate  about  an  orthogonal  set 
of  axes  that  leads  the  blade  azimuthal  angle  by  some  angle, 
<f>,  in  an  analogous  manner  to  the  practise  on  mechanical 
swashplates.   Then,  as  shown  in  Figure  3,  the  apparent 
longitudinal  and  lateral  cyclic  control  matrices  become, 
by  a  coordinate  rotation,  as  follows: 

B       (1)       =      C       (2)    cos    4>    -    C    (3)    sin    <J> 
B       (2)       =      C       (2)    sin    <f>    +   C    (3)    cos    <J> 

where  B(l)  and  B(2)  are  the  longitudinal  and  lateral  cyclic 
control  matrices  respectively.   No  attempt  was  made  to  relate 
these  control  matrices  to  actual  control  stick  motions, 
although  reasonable  estimates  could  have  been  made. 

The  criteria  employed  in  selecting  the  control  lead 
angle  was  the  following: 
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FIGURE  3.   VISUALIZATION  OF  THE  PLENUM 
PRESSURE  LEAD  ANGLE 
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o   Longitudinal  cyclic   . .  produced  negligible 

rolling  moment. 

o   Lateral  cyclic       . .  produced  negligible 

pitching  moment. 


Although  absolute  satisfaction  of  these  constraints 
concurrently  with  one  choice  of  lead  angle,  <f>,  was  not 
physically  realizable,  it  was  remarkable  that  a  single 
value  of  lead  angle,  <J>  =  40  degrees,  provided  a  satisfactory 
solution  in  an  engineering  sense. 

Table  III  lists  the  variation  of  pneumatic  lead  angle 
versus  airspeed  for  satisfaction  of  the  longitudinal  and 
lateral  cyclic  constraints  respectively.   The  selection  of 
forty  (40)  degrees  as  an  engineering  answer  is  in  accord 
with  independent  analysis  done  by  Kaman  Aerospace.   The 
tabulation  of  the  B(l)  and  B(2)  control  matrices  for  a 
pneumatic  lead  angle  of  forty  (40)  degrees  are  presented  in 
Table  IV.   Inspection  of  the  fifth  row  of  B(l)  and  third 
row  of  B(2)  provides  an  indication  of  the  reasonableness 
of  the  solution. 

The  effects  of  plus  and  minus  five  (5)  degree  changes 
in  pneumatic  lead  angle  will  be  described  during  the  analy- 
sis.  This  type  of  sensitivity  analysis  will  provide  an 
indication  of  the  airframe  stability  root  sensitivity  in 
the  compensated  mode. 

C.   THE  FEEDBACK  LAW 

The  study  of  the  impact  of  various  feedback  control 
laws  for  reducing  or  removing  the  longitudinal  instability 
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TABLE  III 
Calculations  of  the  Pneumatic  Lead  Angle  vs.  Airspeed 

B(l)   =  C(2)  cos  <J>   -   C(3)  sin  <J> 

C(2)L 

bi  =    arctan  cirrr 

Li 

B(2)      =      C(2)    sin   <j>      -      C(3)    cos    $ 

"C(3)N 
*B2   =      arCtan      C(2)N 


HOVER  <J>B1   =    42.32° 

»B2    =    43°50° 


35    KNOTS  4>B1   =    37.16° 

*B2   =    39.36° 


72   KNOTS  <j>B1   =    40.87° 

*B2   =    33.09° 


110    KNOTS  cf)Bl   =    39.83° 

*B2   =    36.63° 


130    KNOTS  <J>B1   =    40.61° 

*B2    =    42.56° 
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TABLE  IV 


Coefficients  of  the  Control  Matrix  B 

and  B0   for  pneumodynamic  lead  angle 

of  (40)  degrees 


lx 


HOVER 

35  KTS 

72  KTS 

110  KTS 

130  KTS 

Bll 

-12.96 

-11.13 

-8.474 

-7.593 

-5.121 

B21 

-8.697 

-13.08 

-9.193 

-22.13 

-27.61 

B31 

9.356 

7.565 

6.911 

7.223 

6.726 

B41 

0.0 

0.0 

0.0 

0.0 

0.0 

B51 

1.164 

-1.124 

0.295 

-0.053 

0.174 

B61 

-0.267 

-0.055 

-1.446 

-0.905 

-1.098 

B71 

4.897 

2.816 

2.842 

3.049 

2.030 

B81 

0.0 

0.0 

0.0 

0.0 

0.0 

B12 

5.486 

3.860 

2.086 

0.639 

— — — — — — 

1.711 

B22 

-0.765 

5.042 

2.296 

11.21 

11.53 

B32 

-0.572 

0.084 

0.837 

0.425 

0.304 

B42 

0.0 

0.0 

0.0 

0.0 

0.0 

"B52 

28.72 

22.68 

19.52 

18.30 

16.52 

B62 

-0.106 

-0.643 

0.521 

0.965 

0.585 

B72 

12.94 

10.36 

9.244 

10.60 

9.996 

B82 



0.0 

0.0 

0.0 

0.0 

0.0 
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began  with  an  attempt  to  vary  the  available  longitudinal 
feedback  gains  and  to  investigate  the  eigenvalues  and 
eigenvectors  as  these  gains  were  varied. 

The  modified  plant  matrix,  A' ,  was  developed  in  the 
traditional  matrix  manner  in  the  uncoupled  and  fully- 
coupled  state  variable  format  where: 

A'   =   A  -  Bk 

When  only  longitudinal  cyclic  control  is  considered,  the 
control  effectiveness  matrix  [B]  becomes  an  eight-by-one 
matrix  while  the  feedback  gain  coefficient  matrix  [k] 
becomes  a  one-by-eight  matrix.   The  matrix  product,  B  k  , 
is  an  eight-by-eight  matrix. 


BT   =   [U  ,  W  ,  Q  ,  6  ,  P  ,  R  ,  V  ,  <f>] 


u    w   q    8     p     r    v     4> 


The  earlier  confirmation  that  the  unstable  root  was  primarily 
associated  with  the  longitudinal  airframe  modes  was  the 
reason  for  only  employing  feedback  in  the  longitudinal 
cyclic  control  when  developing  the  modified  plant  matrix, 
[A'].   An  arbitrary  set  of  moderately  damped  oscillatory 
stable  roots  were  selected: 


A    =  1.68  ±  i  8.21 
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which  corresponded  to : 

co  =  undamped  natural  frequency  =  8.38  sec 
C   =  dimensionless  damping  ratio  =  0.2 

Since  it  had  been  established  in  the  basic  plant  matrix 

that  the  short  period  mode  was  the  dominant  instability, 

the  arbitrarily  selected  second-order  system  (a  form  of 

modal  control)  was  applied  to  the  augmented  two-by-two 

matrix,  A' ,  and  the  closed  form  solution  calculated  to  yield 

values  of  k  and  k  .   Gain  values  of  k  =0.05  and  k  =0.22 
w      q  w  q 

were  determined  to  yield  the  desired  results  but  the  applica- 
tion of  these  gains  alone,  to  the  uncoupled  four-by-four 
matrix  and  the  fully-coupled  eight-by-eight  state  variable 
problem  failed  to  produce  favorable  results  and  the  insta- 
bility remained  with  the  aircraft  matrix  at  all  of  the 
calculated  speeds. 

The  search  for  the  acceptable  feedback  law  using  gains 

of  k   and  k   continued  using  the  HP.-9830  with  a  further 
w      q  3 

modified  BASMAT  program  that  would  automatically  search  for 
acceptable  feedback  gains  in  the  range  of  k   and  k   equal 
to  minus  one  (-1)  to  plus  one  (+1.0) .   Although  values  of 
k   and  k   could  be  found  that  would  drive  the  augmented 
matrix  stable  (negative  real  parts  of  the  eigenvalues) 
at  each  calculated  speed,  these  values  were  not  sequentially 
related  to  speed  and,  furthermore,  were  random  in  nature, 
often  changing  signs  more  than  once  as  speed  increased  from 
hover. 


The  decision  was  then  made  to  employ  pitch  attitude 
and  pitch  rate  (kQ  and  k  )  feedback  respectively  based  on 
the  knowledge  that  both  pitch  and  pitch  rate  information 
were  presently  available  in  the  H-2  aircraft.   Programs 
existed  for  both  the  IBM-360  and  the  HP-9  830  for  accom- 
plishing this  search  for  acceptable  feedback  gains  of  k 
and  kQ/  but  while  the  IBM-360  was  much  faster  in  the  actual 
computations  (approximately  twenty  seconds  of  CPU  time  were 
required  for  the  computations  of  the  eigenvalues  and  eigen- 
vectors for  one  set  of  plant,  control  and  feedback  gain 
coefficients  vice  twenty-five  minutes  for  the  HP-9830) , 
the  HP-9  830  allowed  for  a  much  more  convenient  search.   The 
HP-9830  allowed  the  programmer  to  make  "in-line"  decisions 
on  changes  in  the  feedback  gains  based  on  the  previous 
results  with  approximately  thirty  minutes  between  output 
results.   This  removed  the  problem  of  the  long  delays 
encountered  because  of  the  turnaround  time  of  the  IBM-360. 
Typical  turnaround  times  were  two  to  five  hours  depending 
upon  the  computer  usage  at  the  time  of  program  input.   This 
long  turnaround  time  was  a  result  of  the  low  job  priority 
assigned  the  program  by  the  computer  center,  a.  result  of  the 
complexity  of  the  program  and  the  large  amount  of  core 
memory  required. 

The  rationale  of  feeding  back  k   and  kfi  proved  fruitful 
for  values  of: 


k   =  0.45    and   kQ  =  0.85 
q  0 
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or  the  feedback  matrix  taking  the  form  of: 

k   =   [0,  0,  0.45,  0.85,  0,  0,  0,  0] 

producing  stable  eigenvalues  at  all  speeds  calculated  with 
the  exception  of  hover.   These  gains  did  leave  an  unstable 
oscillatory  root  in  hover,  with  a  time  to  double  amplitude 
of  18.15  seconds,  over  six  times  the  three  second  minimum 
time  to  double  amplitude  allowed. 

D.   PNEUMATIC  LEAD  ANGLE  SENSITIVITY  CHECK 

Concern  had  been  shown  over  the  possibility  that 
variations  in  the  pneumatic  lead  angle  (<£)  blowing  could 
cause  dramatic  changes  in  the  stability  characteristics  of 
the  aircraft.   It  was  recognized  that  changes  in  pneumatic 
lead  angle  would  result  in  direct  changes  to  the  control 
matrix,  B.   In  order  to  study  the  sensitivity  to  changes 
in  lead  angle,  new  values  of  the  control  matrix  were  com- 
puted for  pneumatic  lead  angles  of  thirty-five  (35)  and 
forty-five  (45)  degrees.   The  results  of  these  calculations 
are  listed  in  Table  VI.   The  IBM-360  program  was  further 
modified  to  compute  the  new  eigenvalues  and  eigenvectors 
using  the  basic  plant  matrix,  A,  the  computed  feedback  gains, 
k,  given  above,  and  allowing  the  longitudinal  pneumatic 
control  matrix,  B,  to  vary  for  values  of  §   equal  to 
thirty-five  and  forty-five  degrees.   The  computed  results 
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TABLE  V 

Eigenvalues  of  the  augmented  matrix,  A* ,  at 
velocities  of:  Hover,  35,  72,  110  and  130  Knots 


Hover: 


X 

= 

-0.00019 
-0.24473 
-0.47775 
-0.97459 

= 

0.03819 

-r 

i 

0.27676 

= 

-6.78748 

+ 

i 

3.60274 

35  Knots: 

= 

-0.01709 
-0.04164 

= 

-0.24762 

+ 

i 

0.81412 

= 

-0.67857 

+ 

i 

0.13139 

= 

-5.67357 

+ 

i 

3.42617 

72  Knots: 

= 

-0.01445 
-0.07094 

= 

-0.63074 

+_ 

i 

0.37790 

= 

-0.43457 

+_ 

i 

1.47940 

= 

-5.45309 

±_ 

i 

3.07803 

110  Knots: 

= 

-0.00326 
-0.04764 

= 

-0.49037 

+_ 

i 

0.55734 

= 

-0.54836 

+_ 

i 

2.06899 

= 

-5.61868 

+_ 

i 

3.46986 

130  Knots: 

= 

-0.01401 
-0.07035 

= 

-0.39270 

+_ 

i 

0.61107 

= 

-0.59046 

+_ 

i 

2.34511 

= 

-5.31913 

+ 

i 

3.22033 

37 


CO 

m 


o  o  o 

•        •        • 

o  o  o 


o  o  o 

•        •        • 

o  o  o 


o  o  o 

•        •        • 

o  o  o 


o  o  o 

•       •       • 

o  o  o 


o  o  o 

•        •        • 

o  o  o 


CQ 


vo  r* 
o  en  m 

OCOh 

•        •        * 

vo  ^  co 


CTv  VO  CO 
O  rH  O 

r^  oo  en 

•    •    • 

CO  CM  rH 


r*  cm  in 

n  ^  oj 

VO  00  O 

•            •  • 

CO  CM  CM 


VO   "<3«   rH 

•      •      ■ 
CO  CO  CM 


CO  O  rH 

en  co  in 

CO  O  rH 

•        •        • 
CM  CM  rH 


in  r* 

VO 

in  co 

in 

rH 

r^ 

CO 

r- 

CM 

CO 

rH 

in  o 

rH 

r~  vo 

in 

co  *r 

•^ 

r- 

co 

en 

VO 

ro 

00 

r~ 

cm  r~ 

VO 

CM  CM 

CM 

(X)  vo 

VO 

rH 

CM 

CM 

rH 

CM 

CM 

o 

rH  rH 

pq 

•    • 

• 

•    • 

• 

• 

• 

• 

• 

• 

• 

• 

•    • 

o  o 
1  1 

o 
1 

o  o 
1  1 

o 

1 

■H 
I 

rH 
1 

rH 
1 

rH 
1 

rH 

1 

rH 
1 

rH 

I 

rH  rH 
1   1 

CM  «* 

<5r 

r>  ^ 

r~ 

in 

in 

r^ 

CM 

ro 

CO 

CO 

rr  VO 

rH 

>X)  VO 

^r 

m  cm 

C\ 

<Tv 

en 

o 

^ 

in 

T 

rH 

rH  VO 

in 

vo  rH 

CO 

CO  rH 

o 

CTv 

CM 

^r 

in 

o 

vo 

VO 

r-~  cm 

CO 

•  • 

• 

•      • 

• 

• 

• 

• 

• 

* 

• 

• 

•  • 

CO  rH 

rH 
1 

O  rH 
1 

CO 

1 

rH 

O 

rH 
I 

<H 

o 

i 

rH 
1 

rH 

O  rH 
1 

> 

w 
rH- 
pq 

< 

Eh 


CO 


CQ 


o  o  o 

•        •        • 

o  o  o 


o  o  o 

•        •        • 

o  o  o 


o  o  o 

•        •        • 

o  o  o 


o  o  o 

•        •        • 

o  o  o 


o  o  o 

•        •        • 

o  o  o 


o  vo  o 

co  in  <j\ 

r>  rH  CM 

CM  CO  CO 

co  vo  rr 

co  in  r- 

•sp   VO  CM 

in  rH  rH 

co  cm  in 

CO  CM  r-« 

r»  co  co 

m  m  in 

en  en  co 

CM  CM  rH 

r*  r-  vo 

•  •  • 

•  •  • 

•  •  • 

•     •     • 

•    •    • 

en  en  en 

r>  r-  r^ 

VO  VO  vo 

r-  r-  r» 

vo  vo  vo 

CM 


m 


to 


o 

r> 

r- 

o 

CO 

in 

r-» 

rH 

CM 

rH 

CM 

■^ 

rH 

rH 

rH 

CO 

en 

en 

in 

r- 

rH 

en 

in 

rH 

r- 

CO 

CM 

O 

rH 

rH 

r-» 

vo 

in 

en 

r- 

CM 

in 

rH 

VO 

O 

rH 

O 

m 

vo 

in 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

00 

CO 

00 

en 

vo 

o 

rH 

CM 

CM 

rH 

CM 

co 

VO 

r-» 

CO 

1 

1 

I 

I 

• 

rH 
1 

rH 
1 

rH 

1 

rH 

1 

CM 
1 

CM 
1 

CM 

1 

CM 

1 

CM 

« 

CM 
1 

00 

in 

««* 

o 

en 

CO 

o 

-tf 

"3" 

00 

CO 

O 

CO 

rH 

rH 

CM 

in 

CO 

in 

CM 

CM 

vo 

r^ 

CM 

o 

c* 

CM 

in 

CM 

IT) 

-=i< 

Cn 

CO 

r-  rH  «* 

CM 

■^ 

vo 

in 

in 

vo 

.  en 

rH 

CM 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

CM 

CM 

CO 

O  rH  rH 

CO 

00 

00 

r-» 

c-» 

r» 

•^ 

m 

in 

rH 
1 

rH  rH 
1   1 

rH  rH 
1   1 

rH 

1 

1 

I 

1 

J 

l 

i 

» 

I 

1 

en 

en 

en 

en 

en 

Eh 

Eh 

Eh 

^ 

w 

o 

0 

o 

W 

o 

0 

o 

tt 

0 

0 

o 

« 

0 

0 

0 

« 

o 

o 

0 

> 

in 

O 

in 

in 

O 

m 

in 

o 

in 

in 

o 

m 

in 

O 

in 

o 

CO 

*tf 

^ 

in 

CO 

*tf 

■^ 

CM 

CO 

"^ 

•^r 

o 

CO 

«* 

** 

o 

CO 

^r 

«* 

K 

CO 

r^ 

rH 

CO 

38 


showed  that  the  variation  of  the  pneumatic  lead  angle  by 
plus  or  minus  five  (5)  degrees  had  very  little  effect  on 
the  stability  characteristics  of  the  augmented  plant 
matrix. 
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IV.   CONCLUSIONS 

A  study  has  been  made  of  the  basic  stability  traits  of 
the  Kaman  Aerospace  Corporation  XH-2/CCR  helicopter,  which 
is  presently  being  constructed  under  NAVAIR  contract  as  a 
technology  demonstrator  for  the  Circulation  Control  Rotor 
concept. 

The  airframe  was  defined  by  contractor  generated 
stability  and  control  derivatives  which  were  then  used  to 
develop  eigenvalues  and  eigenvectors  for  the  system.   The 
plant  matrix  (which  characterizes  the  airframe)  was  generated 
by  the  MOSTAB  program  as  modified  to  accommodate  the  CCR 
system,  and  the  matrix  coefficients  represent  the  output 
from  the  program  when  it  was  operating  in  the  18  degree-of- 
freedom  situation,  i.e.,  six  airframe  degrees-of-freedom 
plus  flapping,  torsion  and  lead/lag  degrees-of-freedom  for 
the  four  blades  (6+3x4  =  18).   The  uncompensated  airframe 
characteristic  roots  at  hover  are  in  close  accord  with 
results  obtained  by  the  contractor  for  the  six  degree-of- 
freedom  airframe,  including  a  mild  non-oscillatory  insta- 
bility (t2  =  37.1  sec.)  that  has  been  identified  as  due  to  a 
longitudinal  short  period  type  mode. 

A  feedback  law  has  been  identified  using  pitch  attitude 
and  pitch  rate  feedback  into  the  longitudinal  cyclic  control 
that  provides  reasonable  characteristic  roots  for  the 
airframe.   Presumably,  further  improvements  could  be  obtained 
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by  providing  feedback  in  the  lateral  cyclic  control  system. 
The  effect  of  varying  the  cyclic  control  lead  angle  on  the 
pneumatic  swashplate  was  investigated  and  found  to  be  slight. 

It  should  be  noted  that  a  unique  feedback  control  law 
is  not  possible  in  modal  control  theory,  when  multiple 
control  inputs  (longitudinal  and  lateral  cyclic)  are 
available.   Another  way  of  stating  this  fact  is  that  it  is 
possible  with  several  control  inputs  (and  feedback  laws)  to 
have  the  same  eigenvalues,  but  with  different  eigenvectors. 

Finally,  the  characterization  of  the  eigenvectors  and 
identification  of  the  eigenvalues  with  relevant  modes  was 
made  possible  by  using  a  form  of  root  locus  analysis  where 
the  prime  parameter  was  the  amount  of  cross  coupling.   The 
trajectory  of  the  characteristic  roots  as  cross-coupling 
was  linearly  altered  provided  a  physical  insight  into  the 
history  of  the  various  roots. 

Future  studies  are  suggested  to  include  estimating  the 
mechanical  gearing  to  the  cockpit  controls  and  then  obtaining 
airframe  response  time  histories  for  selected  control  inputs 
such  as  stick  doublet  type  motions.   Time  histories  can  be 
generated  quite  readily  using  principles  from  control  theory 
combined  with  calculated  relations  for  system's  transition 
matrices.   It  is  quite  possible  that  time  history  calcula- 
tions for  the  compensated  airframe  will  provide  a  better 
guide  for  selecting  the  control  law.   In  any  event,  such 
studies  will  enhance  satisfaction  as  to  the  question  of 
airframe  response  behavior  being  reasonable. 
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0C7Q  2023     F"IPM£(      (  11-0,  5X,31Hr'-»P    MOdFTFO    PL->N"    'M-3^,     i,-3G,/) 

0030  EviC 


43 


0001 


SUa9?'rr  Nc    CHSFQ(4,"I.  C,*I9*) 

THIS     ;i.'m-ii)TT  \p     FINOS    THE    CFPF  i  r  Lt-MTS     OF    "*  t-c    r  MS  '  ACTE~  I  ST  T  r. 

P"CVNV<IAl'    U5IMG"TH-    LEVEIES    iL^""'.!TMM i_ 

Cnw-^CN    ZED  (10  ,10.10) 

oimf^v  in  .\(  is.LJi .:  u  n  ,:;e«p  i  is  ,10) .  go  en  a  o,  101 


0002 
3003 


0004 
0005 


DATA    AT?mo/  100»0.0/ 
100O    FH7M4T    (  1H0.5X,31H*HE     «",  T  5 1  x    CrEFFTTI ENTS    OF    THE 


3  3H.1|tJMco  i  rr 


=h€- 


iLVF'l1-    '«1~»IX 


0006 
3  00  7 


I  00  I 


F1R-4  4T     (IHO,        5X.29HTHE    MATRIX    COEFFICIENT    OF     S**,T1/1 
FOpvar     (JJliE_2_Cj_Z_L 


L003    nvu*     (  LFC.  7i5(  1H*)  > 

CALL    CHQcri (i  ,N,C) 
01     65     T -L  ,\ 


000  8 
00O9 
0  313 


0011 
0012 
0313 
0014 

001^ 


r-5 

70 


A' £"•">(  I.  I  )  =  I  .0 

oi  3  0  :  =  i  ,n 


DO     'SO    .l=l,f 

30     ZFCtN,  I  .  J  l=lT="D(',.i  ) 
TF     (N3".NK.3)     Q2     "     7  1 


UO  TT  E    UTf  IC03  ) 
W='~E    (6,1000) 

."!='.- 1 


0016 
0017 
331P 


0019 
0020 


3021 
0022 

0023 


W3.  IT?    (6  ,  100  L  )    ^ 
0~    3«5    !=1  ,\ 


55   w°  r 


( 6,  I 


i — 1  iTCMo  (  ;  ,  j>  ,  j-l  ,»o 


UTT+ 

002  5 

3  326 


71 


C"1    tO    1  =  1,  N 

0  "■    40     J  =  1  .  - 


c~   10   i=i, r> 

MMMsr,|-T 


002  7 

0028 

0030 
QQ3  1 


TF     (*.£■}.  I  )    GO    Tl    55 
IF     (M'M.NF.O)     C-C    T"1    60 


^w<^ 


I  0  .  1  q  3  I  )     VIMN 


0032 
0033 
3334 


C"    45    J=1,N 

WP'fF     (^,1002)      (  A,TC-<P(  J,K)  ,*=!  ,N) 


o°  9  3  i  :  =  i  ,\ 

i 


37T 


-a 


'3    ,J=1 


0035 
0036 


0037 
0033 
0039 


90 


z=c?np, T i, j )  =  atfmp{  :  •  , J) 

0"  15  J=l,': 


O'l  15  X-l. 
PRG3I  ),<) =Q.3 

nn  15  L=  1.: 


004  0 
3341 
0042 


«:i5(J,K)=c5':oU.iO*(  AU.L)*&TEMO{i_,i<  J  J 
D"1  13  J=1,N 

D1  13  K: 


TT 


-L- 


0143 
3344 


0045 
0046 
334-7 

Onoi 
3332 


13 

55 

-to- 


ATJMOt j,k ) =F3CC( J,<) 
Dn  13  J»-l'.N 


?TE«P(J,J)^~g*fMJ,J)»C(>i-!»l> 
RFTtj3N 

ENC  


0003 
0004 
-W-3-5- 


C006 
0007 


0008 
0C09 

0010 


001  1 
0012 
331  3 


0014 
0015 
3316 
001  7 
0013 


00L9 
0020 
0021 


0022 
002  3 


3  32  4 
0325 


0027 
0028 
0029 


3030 

oom 


3  33  2 
0033 
C03* 

3335 
0036 


SL'PSHjTTNE  CHFECACi,  N,C  ) 

h-^-nc  i r -.m  j(  i  i)  .C  (  11  )  ,fi(  13,  10)  ,M  10,  1Q)  ,0(300) 


20  C  t  I  »  -0.0 


r  r  n  n  )  =  i .  o 

nn    14  M= i , ^ 


K  =  ■" 

L  =  l 


J  (  L  )  =  J  ( l  J  *  ! 

i  L  -  * )  3 . ; 


J 

O'l 


TFT 

VI" 


s-f-T-t 

n  )=j(  u+i 

10  ' = 1 . v 


10 


B( 

K  = 


0  *K=  l,« 
=  J<  T  ) 

■J  (  k-  «  ) 


I.Kk  )  =  i( MP ,NC ) 

KM 


p  (  m  )  -Oc~  (  R  ,'4  ) 

CO  6  '  =  l  ." 

T=  ( 3TT  >-('!-" +L  )  I  1,  i,  5c 
ton-  t  nijf 

VI  sN-V-H 


14 


I  '♦  I  =  1 .  K 
M )=C( M  I +C ( I  )*  (-1  .0) 


P"TsT  2300 

F'  -viT  (  1HC  ,5^  ,1SH  =  - 

-  E  "iu"i 


IN  CHRFQA) 


44 


oooi 


TjTTNf= 


~rr-r 


i*'G 


s T  :n  :  (  io  f  10 ) .  3  ( 10 1  i  o  j  ,  xo-»~  (  i  i  > ,  i  tn  v  (  ic,io).*(ii) 


C0G2 


-4U4- 


000^ 
'005- 


000"" 
O^OJ 


-  =  I 

1      !  =  I   r  K  C 


0010 
001  1 


I     S( 


6  T 


V/<  '  .  J  )=C 

'    ,     ))    ='    (  T     ,    J) 


0012 

001  3 


0014 

0015 
001* 


2   :=l 

!W(  i .  I  )  =  1 

sun 


U4 


t-U- 


NO=0 


-r-=- 

?  F 


a"    ( Q ( y  . : )  )--is ;   ( crM»  i )  r,5  .*. 

=  c  [K    I  '    ' 


00S2 


45 


000  1 


SJS^n'jTi  ,\j=     5T-4ST    (N,  \  ,Eio»  .  E  r  ■';  I  .TK  HW) 
rml<;     «U3C1UT!NF    ncTE°'-l  i  N  FS     ThF    5T4.7F    "3ANStTt0M    M4~RIX    USING 


g,-?<  lo >  ,ft",t  do)  .  s?s  in.  13  ) 


0002 
0003 


Cj"M.1N    CHI  (  10.  10,  10) 
I'VJENSCIM     MIO.IO; 


0004 

0035 
0006 


cn«oLcX  r.i  <  io.  10)  ,cak  io,i  ji  ,Ci2(  n,  lot  .TC4I  io,io>  , 

*       CENOM(  10)  ,CETG(lOf 


1000    F<!**AT     i  1-*0.[3X  .25"~wc     g[r''CAlTS     IP    ~HE 


mai 


13H    ~5i,'.|?TT'0N    MATRIX 

FOOM  i'llHQ,         ^X,25HTMp     MA 


-**-e- 


) 


f:?»mt(  IP8F.L5.4) 
PQa.<lA~(  1MQ  , 5X,25H~HF.    MdTOIX     r.lFFFICTFMT    OF 


0007 
0008 


1002 
1-003 


000*5 


*       5M    EX?< ,  l°E  12.6, 7H I T*S'M(  , i?c 1 3 . 6 , 2H ) T  /  ) 
1004    FTSMAT  (IHO,       SX.25HTHE    '•i.iT^lX    COEFFICIENT    TF 


5M  c  XP<  .1PE13  .ft  ,    gtttr-/ ) 


0010 

T3T2 

0P13 
PC  14 


1005 


ft3m4t   (  1hc,45(  ih*)  ) 

tc  (  :kmow.nf.o  )   go  t^   300 


PS  TMT    100? 

300     Pr    10    K=l,N 

c^ir,  (k)  scv'flx  (=' 


IKU  )  .  FTGI  (K  )  ) 


0015 

0016 

0017 

0013 

CC1Q 

TJTT 

002  1 

0022 


10 


D°     10    L=1,N 

Ci(K.L)     =    CuPLX(A(K,L).O.OJ 


«-^ri 

T«={  IKNvOW.NP  .C)    GO    r"     7C0 

pg  t-;t  loop- 

70  5  on  i;  k=i,\ 

15    RP\-'MK  )=CETG<  T  )-CE!G(«  ) 

pa   500   j  =  i  .v 


0023 
0024 
0025 
0026 
002  7 


100 

300 
400 
TTT 


IF    (  J-T  >     ICC5CC200 
TF     tJ-1)     I  10.  1  10, i 50 


IF      I 


■1)     3^ 


*39- 


T3TT 

0029 

0030. 


<tjb.i 


TF    U-I-l)     11C,  110,  150 
TP     (J-r-l)     110,150,130 


DTI    5    «  =  l  ,v! 
rn    5    L=l.N 


5    C  i  11  K , L  ) 


■U- 


0031 
0032 
0033 
0034 
0035 


on    2  0    K  =  1 ,  i\ 

Ca  L(K,K)=CMK.K)-CE!GU  ) 


D1-1    20    l-I.N 

20    CMIK.L  )=C^1(K,L)/0E,J0M(  J) 

go  tc  5oo 


0-  40  K=1,N 
C?  40  L-l.N 
C'  2(  K.L)  =C  M  K 


0036 
0037 
-0038 


150 
-40- 


<K  .1) 


•003Q 
0040 
004  1 
0042 
0043 


DO     25    K  =  i.N 
C\2(K,K)=CMK,«)-CE!3U) 


0°    2' 


Tv*- 


3  544 
0045 
DQ46 


25 


CA2(K,L)*C* 2(Ktl >/DENQM( J) 
C"1    30    K=l  .* 


C"i    3  3    L  =  l.»\ 

TGMK,L)  =  <  C  .O.C.O) 
0°    30    H=l  ,f\ 


0047 
0048 
0C49 
0050 
005  1 
TTOTT 
0053 


30 


co 


^(  K.D  =TCA(  K.L)     «■    Ci  K  K,V)*C-'2<  *.l  • 
35    K=L,f. 


»q  l  »  i  h< 

KK.LI  ='Ci  (K  ,L  1 
MTNUF 


35 

50  0 


r  a 


"TT- 


(il'lir,  (CSri'Jd  )  )  I     4^,:>J,4b 


3055 
0056 
00r*7 
0053 
0  05  9 

rrTTFn 

006  1 
306  2 


-*ttr 


P3 

TTT 
XH 


K"lf?w.NE.3)    GO    ro    83  1 
T     100  1  .     ETC3(  I  ),  EIGI  (  '  ) 


5    1  =  1,*' 

K.L) =    =C*L  ICA1 (K.L)  )*2.0 


5  ( 


6    K=l  ,h 

6L=  l,N 

IV.,i  .CUSPS  (K.L) 


u 


0063 
C064 
C06  5 
3366 
0067 
TFTFT 
0069 
0070 


66 


rn 

T  P 

DC 

po 
no 

-S-2. 


TMUF 
KNTW.')5»C)    GO 


8C2 


T    10C2  ,     ( S»5(K,L J  ,L=1  ,M> 
T    1003,     ElCPJI)tFIGI(I) 

5    K  =  l  ^ 

5    L  =  l,N 
K.i.)aAI'<:CtCll(ft.L))*2.3 


30 
55 


0071 


OT    5  6    K= 1 , S 


46 


0072 


DO    56    L=l^ 

C  H  U  7 tKiU^^MK.I.  ) 

cIntTnue 

»p  {'.  KNGW.NP.J)    GO    Tn    600 


0074 
007  6 


007  7 


0373 
0079 


35 
*=*    ' 
210 

770 


pu 


CUL4T=    CC,J,t»L  =  X    6!3ENVEC"}5     («    MT0!  F  '  CA  t  ;  7^  ) 


o^ f  -it — loiv 


DC    221     K=1,N 

r->t  *  .  1  )=SC3'C"HI  (I  m.k  .  1  )--*?     »    CtjT  (  T  .<.  I  )*»?  ) 
?=<;  1  K.2  )  =-4  7  JN2(CHt  (  !  .K,  I  )  ,  CMt  I  IM  ,  K  ,  1  )  )     »(  -5V.2 
P:;f!*      iai7.?P<(K.II.<PslK.7l 


51  NT    l017,SP;<K,L),cP^tK.2) 
•"■i-T"'JP 

C5Md7(  1H0.5X  ,26HrHP     COMPLEX     EIGENVECTOR    IS  ,  /  3  X  ,  9HM  AGN  I  T'JOE  , 

13X,11H?HASE     <OE~,)/l 
'^"iTI  L  J  ?!=2C  .7) 


008  L 
0082 
0083 


-3*4- 


0084 


003*^ 
Q0S6 


1006 


I 


100  7 
:    **    ENQ 


F    M--0!  P!  :iTTCN. 


50     IF(  I  KNCiW  7NE  .0  )    GG    ? 
oi\\r    1004,    r.'G=(!) 
-A44 — P^3 — 6  0     K  -  1  .  N 


0037 
0038 
0039 


304 


0O90 
0091 
009  2- 

C09"3 
009^ 


60     S 


DO 

CM 


6  0    L=l.l^ 
SJKfUagSAL     (C^IU.D) 


ft  I    «-  1 .  f 
61    L=l  »i* 


0095 
C"»9* 
0JC7 


61     C"iNTrilJp 

IFtrKNnw.NE.0)    GO    TH    600 
pq    7  5    y-ri,\ 


*NT     1002  ,     ISP"  (K  .L  )  ,1=1  ,VJ  ) 
(I.GF.iXJ     =>CTL*"-l 


C093 
0099 
0  10  0 
0101 


75 
600 


P^ 
IF 

T= 

CM 


7C    700 


47 


0001 


=)0  0  UT  f  N  E  P  ^  r  C  T  (  N  ,  \  ,  ;j  ,  V  ,  T  }  ) 

>iip.?-!  ittNc    :;cc;    s    ^nnrFTPn    ELS^LSIOU   """^    t^    EJHQ   ~hf 


0002 
0003 


TS    OF  "d    Ffi.YN  "3MI  ST  . 

*.=  NS'1N    4(20)  ,U(231  »V(23).H(21)  .3 (21  )  ,  C ( 21  ) 


0! 


0004 
0035 
0006 
0007 
3U3P 


NC 
00 

I  HI 
P- 

Q= 


11=1,  MC 
I  >  =  at ? ) 
c. 

L- 


000<? 
0010 
OOl  I 
0012 
0313 
00 IV 
0015 
3  Tift 


D  = 

3    TF 


(H(  1. 


) 4,2,4 


V( 

U( 

TT 

100?    H( 


NO- 
NO  = 
Trtnr 
:  i  =h 
O  3 


0. 
3. 


n  +1 ) 


001' 

OOlfl 

001^ 


4  T  = 

5  IF 

6  3- 


(NC 
(NC- 


5,  13 
7,6, 
W-54 


0  . 
7 


0020 
3321 
002'2 
0023 
0  324 


7  TP 

-* — B^sr 
Q  = 


50 
(MC- 


5,6.' 

t  J ) 


H(  1)/M(3) 

Ti  70 


00  2  5 
0026 
0027 


q 

13 


r  p 


(T3S 

6V  = 


R  NC-l  )/HMC  )  )  -i  a  5  (*<  2  >/M  1  )  )  )  10,  10, 19 
FV 


G02* 
332q 
0030 
00  31 
J  33  2 


oc 

ML 

f- 

h( 
JJ HJ. 


11  !  = 

a  I  C  * 

mr 


NL)=H{ T ) 

I  1  =P 


0033 
0O34 
W3-S- 


12 


T  C 
P  = 


p  I  13,  12,1 

H+-5 


0036 
3037 

00  3* 
003O 
3  34  3 


13 


-rr 

16 


0  = 
—re- 


,/0 

rrrrrr 

,C-1Q 


T*" 


Q  = 

c  =  i 


0041 
0042 

0-O4-3- 


-SM- 


NC  ) 
)  »H  (  \r  ) 
»1)=0. 


0044 
3345 
0046 
C047 

3  34* 

004  0 

0060 
0051 


-rrr 


C( 

MP 

-r—r 


NC 


r-i 

j»t ;  ioao 

I  1  =  t   .  N  P 

L 


I  = 


?1 


0052 
3  36  3 

0054 


i ) =h( :  )♦  =  *"( * *i ) 

21    C( T) a8(T  >+b*C(t  +  l ) 

t  c  (  -.  a  <:    (  a  t  1  )  /  u  ( 4-4-4- 


R  = 


(2)  ) 23,22  ,23 
1. 

?0 


50.50,24- 


48 


005* 
3356 


23 
_3_3_ 


B( 

-G4- 


8(1 

Li 


!  i 

i(T  )-P*f(  I>1 )-0*BU+?) 
2)  )  32.31  ,32 


0057 
0058 
-W-5-Q- 


-*X 


=  H( 

^£4 


3363 

0061 
006? 
3  06  3 

006  4 


31 

33 
34 


I 
TF 

T  C 


Z>  )  32.31  ,32 

*  t  8(  2>/M  l))-c  )  33.33,34 
J  (?(?)/>■  (2  i  >-r)  33.  33.3* 
"    ( 3(  H /H(  1)  »-=  1  70,  7 C, 34 


<  Z) 

ss 


0065 
0066 
3367 


r  = 

P  s 


?-r o J -*r (4) 


(  0  I  i  6 ,  3  5  ,'  3  6 
B-i. 


0069 
0069 
3070 
0071 
007? 


"3^ 


0  = 

0" 


'  (0  +  L  .  ) 

!49 

■(  ■  (  g  >  ''C  (  3  )  ^8  (  t,  I*"  I  frl  1  /"* — 

M-B<  2)  >C°:'?*9  (L5*C(3)  I/O 

TiLir. 

H7T. 

)20 

iC=4 


3373 
0074 
0075 


►  9    C 


C  =  ( 

FT! 


S4- 


3376 

0377 

0073 
007* 
0080 
008  L 


V(NC  )=0, 


51 


IF 
U( 

u( 


if1 

1  =  3 


151 

./3 


52,52 


52 
-54- 


■^~ 


3082 
0083 
003  4 
3085 

0036 


54 


H( 

GO 


^1  [+i  ) 


0087 
3  38  8 
0G89 


71     CP 


3=  VI  71  ,72  ,72 

43- 


/(Q*2.C) 
73 


7? 


po- 


0090 
009  1 
COP? 
3393 
0094 
C095 
03<36 
009" 


73 


PP 

U( 

'J' 


/2.0 

D  |  a* 
)=-° 


2-C^ 

7«s  m 

-op 


'4 


tct — mr 

(NCU  J 


V 

V( 
-GO 


+  1  )  = 
)  =  -v 
is 


0099 
3399 
3100 
CL01 
0102 


75 

90 


'M 

L( 


D)  ai 
+  1  )  = 

+  L)  = 

TNUg 


,  3  C  f  8 1 


0  10  3 
3134 

0105 


81 

32 


-(cc/apc    (op))*(49S    (»o)*S3'T    (M) 


vT 

U( 

-U4- 


>=Q? 
1 


/tl(  \C  +  1  > 


0106 
3137 
0103 
0109 
0110 


^6 
77 


!=1  . 
=  ^(  I 

A 


HZ 
«-2> 


133 


3'l 


49 


APPENDIX   C 


MODIFIED    "BASMAT"    PROGRAM  FOR  HP-9  830 


CIM 

■ttr 


A  <  E, 8), 0(8.8) 

Mta.n.t  t  i,a) 


8ASQ0010 

3AS00030 
_8A  SO  0.040 
BA SCO 0  50 
8AS00060 
BASCCC70 


CIM  Z(8,3) 
FI>EC  6 


7  CATA 

8  CATA 

9  CATA 


10  CATA 

11  CATA 

12  CATA' 

13  CATA 

14  CATA 


8ASC0Q80 
8AS00090 
9ASC0100 

BASCOUO 
8ASC0120 

"8A~s~<rcno 

8ASC0140 
3ASC015Q- 


15  DATA 

16  MAT  REAC  0 

17  PRINT  "CCR  130  KNOTS  THE  BASIC  A  MATRIX" 


18  MAT  PRINT  C 

19  PRINT 

20  MAT  ESAC  M 

21  PRINT  "THE  BASIC  CCNTPCL  MATRIX  IS" 

22  MAT  PRINT  M 


8AS00160 
BASCOUO 

8ASC018& 
8AS0C190 
BASC02C0 

~B~ATC0  210 
8ASC0220 

-BASCC2  3-0 


23  FOR     I5=-l    TO    1    STEP    C.  1 

24  FCR    J5=-l    TO    1    STEP    C.l 

-25    YU,l)  =  Y{l,2)  =  Y(l,5)=Y(l.M~YU.7)=Y(lt8)=Q 


27  YU,3)=I5 

28  Yd  ,4)=J5 

29  MAT     £-M»Y 

30  MAT    A=C-Z 

31  G1=C 

"32  GCSUB    2C00 

33  N  =  8 

i-lNT"-IhE    REWSEC    PLANT    MATRIX    AT    130    KTSt     WITH    IC(Q)^«tl5 


BAS00240 
BASC0250 
8ASCC260 
8ASC0270 
BASCC28C 
"BASJT0290 
8ASC03CO 

BAsooaaa- 


35  PR1NT"ANC  MThETA)  =  ",J5 

36  MAT  PRINT  A 

37  RCO  IM  K  (8,0) 


3ASC0320 

8ASC0330 

-t»ASC034O- 


38  PRINT 

43  FOR  1=1  TO  N 


44  F(I,1J  =  M1,1) 
46  NEXT  I 
-46  FCR  1  =  2  TO  N 


8ASC0350 
BAS00360 
~3ATCU370~ 
8ASC0380 
aAS00390 


'.Q    1.5=1-1 

52  FOR  J=l  TC  N 


)4  F(J,  I)  =  0 

c6  FOR  K=l  TC  N 

58  F(J.I)=P(J,1)*A(J,K)*F(K,L5) 


8ASC0400 
8AS004  10 
BA500420- 
BASC0430 
3AS0044G 


~6C  NEXT  K 

61  NEXT  J 

62  NEXT  I 

63  FCR  1=1  TC  N 

64  FCR  J=l  TO  N 

65  K(  I  ,J)»M  I,J)- 

66  A( I  ,J)=F( I  ,j) 

67  NEXT  J 


5A  5W4  50 
8AS00460 
SASCa4  7q 


8ASC0480 

3AS00490 

BASC05-eO- 

BASCC510 

8ASC0520 

"3a  saa  530 

3AS00540 
^Ai-OOiSCL 


"6  8  NEXT  I 

69  GCSU8  1CQ0 

70  IF  £  NCT  EULAL  TC  ZERO  TFEN  76 

"il  GCSL8  2GC0 

72  PRINT  "MATRIX  IS  UNCCNTPCLLABLE" 
74  GCTC  HO 


BASC0560 
BAS00570 
3A5C058O 


76  FCR  1=1  TC  N 
78  FCR  J=l  TO  N 

"ts-qt  n^  r=  c 

80  FCR  K=l  TO  N 

dl  flli  ,w)=fit  I  .J)*F(  I  tKL»P(y,J1 

28  NEXT  K 

32  NEXT  K 

84  NEXT  J 


8ASC0590 
BASC06C0 

"3ASCC6  10 
8ASC0620 

-3ASuQt3-Q 
3AS00640 
3ASCC650 

-4A-5-0G^60- 


86  NE>T  I 

e8  E1=C   

^-c-Fxnr-i  =  i  TO  N 

92  FCR  J=l  TO  N 

93  A(  I  ,J)=K(I  ,J j 


94  IF  11-^J)  NCT  rCUAL  TC  ZSPC  TMEN  100 
96  E2=ABS(C(I,J)-1) 
If  GOTO  10 


8ASC0670 
3ASC0680 

"3TrS"C0690 
3ASC07C0 
3ASC0710 
8ASC072tr 
8ASC0730 

_BASCQ7_40, 
3ASC0750 
3ASC0760 

-3A^S0077O- 
BASC0780 
8ASC0790 
JASGOflCO 
BASC0810 

.3AS00820 
3ASCC830 
3AS00840 

-SA-S£&*50- 
3ASCC860 
8AS00870 

BAseoaeo 


IOC  E2  =  AeS(Ct  I, J)  ) 

101  IF  AeS(El)>AeS(E2)  THEN  104 

102  E-1=-A8S  <E2) 

103  GCTC  1C5 

104  E1=A8S(E1) 

105  NEXT  J 

106  NEXT  I 

Jj0  7  PRINT 

"IOtTF" " (El-lE-CSXO  THEN  Lie 

1C9  PRINT  "PLATN1  IS  NUMERICALLY  UNCONTROLLABLE,  OEV  I  AT  I CN  =  " ; El 
UC  GO-SI  E  24G-G- 


112  PRINT  "OPtN  LCCP  CALCULATIONS" 

114  PRINT 

116  PRINT  "CENCM  CGEF  IN  ASCENDING  PQWgftS  OF  S» 


50 


12C  GCSL8    2199 

122  FCP     111    TC 

124  PKlNTTTTr 

126  1(1  >=cm 

12a  NEXT    I 


N6 


8ASC0890 
3ASCQ900 
SAS00910 
BAS00920 

-aA-sso-9-3e 


liC  GCSC8  3399 
132  GCSUe  5000 
*3C  next  JS 


8AS00940 
BASC0950 
BASC0960 
BAS00970 

JASC0980_ 
8  A  SCO 9- 90" 
BAS01000 

-3A-SC1010 


44  C 
45C 


NEXT 
END 


15 


TTTCC  N7=l 
10C2  £=1 
1CC4  FCP  1  =  1  TC  N 


ICC6  FCR  J=l  TI  N 
10C8  PU,J)=0 
1G10  G(  I  ,J)-»AU  ,J) 


BASC1020 
8AS01030 
8A5C1040  - 
BAS01050 
BASC1060 
3ASC1070 
BASC1080 
-BASCIC-SO- 


1012 
1014 


NEXT 


1016  FOR  1=1  TC  N 
1018  X(I)=1 
-1X-1-9— P-I  I  ♦  I  i  =  1 


TC  N 


102C  NEXT  I 

1C22  FCR  1=1 

1024  C-C 

1026  K=I 

1028  IF  (A8S(G(K , 1) )-aeS(C) )  <=  0  THEN  1034 

10.2C  C=G(K,T) 

1032  N7=K 

-1034  K  =  K  +  1 


3AS01100 
8ASC1110 
8A5C1120 
8AS01130 

BA I  CI  140 
BAS01150 

BASC1160 
8ASC1170 


1036 
1038 


IF 
IF 


(K-N)  <=  0  THEN  1028 
IG(N7,I))=0  THEN  UC2 


104C  IF  (N7-1XC  THEN  1102- 


3AS01180 
3ASC1190 
8AS012QO 
3ASC1210 
3ASC1220 

SAsSIzaS- 

3ASC1240 

_8A-S042-50- 

3ASC1260 

3ASC1270 

3AS01280 

BA5C1290 

BAS01300 

""8"A"S"CT3TO~ 

3ASC1320 

-8AS01330- 


10<2 
1C44 


IF  (N7-I)=0  THEN 
_F_C_R_  M=l  JTC_  N 


IC66 


"1 04~6  T=G(I,MT 
1C48  G(I ,*)=C-(N7,M 
1.CS0  GIN7,^)-T 


1052     T=P(I,fJ 

1054    P(  I  (M)*P<N7,M) 


•test    P(N-7,*)-T 
1C5E    NEXT    M 
1C6C    T=X(I) 


TCT2~TTn='XTTr7-J 
1C64    X(N7>=T 
-1C66    XI  I)=XtI)/£  II, I) 


1066 
1C66 


T=X(I) 
T  =  GU  ,1) 


1070    FCR    M*l     TC    N 


3ASC1340 
8ASC1250 
3ASC1360 
BASC1370 
8ASC1330 

"3A:SG1J903; 
3ASC1400 
3ASC1410 
8ASC1420 
3ASC1430 
BASC1440 
8A~S~Cl450 
8AS01460 

-3ASC1470 


1072 

1C74 

TGT6T 


P(I«M)*P( I ,M)/T 
G<  I.^)=G(Itf )/T 
TTETT — 


~W 


1C7S    FCR    J=l    7C    N 
10E0     IF     (J-I)=0    THEN 


1096 


1082  IF  td(J,I)-C  THEN  1C96 
X  (  J  )  =  X  I 

T=G(J. 
FOR  N7^ 


1C64 

10E6 

TCT6" 

IC9C 

1C92 


(J  )-GU,  I  )*X(I) 
I) 

TP  H" 


FOR  N7=L 

P(J,N7)=P(J,N7)-T*F(I,N7) 
G( J.N7}-C(J,N7)-T*C  II  ,N7) 


1C94  NEXT  N; 

1C96  NEXT  J 

1C5E  NEXT  I 

LlCC  PETLRN 

1 1C2  PRINT  "VA  TRIX  IS  S  INGllAR" 


8ASC1480 
3ASC1490 
8A3015CO 
3AS01510 
3ASC1520 
~5A~S7J  r530~ 
3ASG1540 


11C4 
UC6 
15CC 


15C2 
1504 
15C6 
2000 
2CC4 

3200 
-3201 


¥=TJ 

RETCRN 
-FXR— L- 


1  TC  N6 


(  I  ) 


R (I )=C 

NEXT  I 

RETURN 

PRINT  ••*»»**<*»»»?»»«<«**»««»**««**«■*»*■«****•«**•• 

RETlRN 

N6=N+1 
CI*  Ci 


100) 
,1  TC  fi6- 


3ASC1560 
8AS01570 
BA5C1580 
BASC1590 
8ASC1600 

sasoisto- 

3ASC1620 
BASC1630 


32C3 
32C5 
I2C7 
32C9 
3211 
TZIT 
3215 
4X 


CU>  =  0 
NEXT  I 
C  (N6) 


i 

f=l    TC    N 


3ASC1640 

6ASC1650 

8AS0I66O— 

dASC1670 

3ASC1680 

"BXS~CT690~ 
8ASC17C0 

-3AS0171Q- 


FCR 
A=C 

r=T — 

J(  1  )=l 
GCTC 


2-2-L 


3219    J  IL  )  =  J(L J  +  l 
3221     IF     (L-M»=0    THEN    3225 
■223     If     (I    *)>Q    Tt  CN    3277 
3225    Ml  =  »*-1 
3227    FCR     I  =  L    TC    y  1 
3229     11=1+1 


3ASC1720 
3ASC1730 

-9ASC1740- 
3ASC1750 
8AS01760 

"BA3"C17  70" 


51 


3221 
3223 


1225 
5237 


J(I1)=J(I)+1 

Ji£XJ— i_ 


FCR 

FCR 


2239    NS=J (I ) 


1=1    TC    f 
Kl  =  l     TC    I* 


8ASC1780 
BASC1790 


BASC1800 
BASC1810 

-&A-5^Hr&20 
8ASC1330 
3AS01340 

"BXSXI 8  50 
8ASC1860 

-flASCia7Q 


32*  1 

2243 


Nl=J(Kl) 

e(  I,K1)=A(N5,N1) 


3241>    NbX  I    KX 
2247    NEXT    I 
JL2Ai-JL=JLtL- 


5251    GCSU8    3299 
2253    C  (A  )  =  0 

3255  FCR  1^1  TC  P 


BASC1880 
BAS01890 
-8-4  SCI  9  00-- 


2257 
2259 


L*M-I+1 

If  (J(L)-(N-M»l)  KC  THEN  3219 

rnTT=T7*=TWCTJ  ;  U  t  F  FN  32T7 

I 


CT- TF 
3263  NEXT 
12,65  h2=N-M+l 


3267 

3269 


FCR  1=1  TO  A 

C(y2)=c(f2)+c(i)*(-i)ExP(y) 


32  71  NEXT  I 
32  73  NEXT  M 

3275  RETURN  _ 

J277  PRINT  "  ERRUR  IN  CFR  ETA""" 
2279  RETURN 
-3ZS9  11=  C- 


BASC1910 
BAS01920 

"B~ASCX93C 
BASQ1940 

JULSfliaiQ. 
8ASC1960 
BAS01970 

-3A5C1980 


33C1  FCR  1=1  TC  H 

33C3  K=I 

33C5  IF  ( e ( KT 1 J  )  NCT  EtUAL  TC  ZERO  THEN  3313 


8ASC1990 
BASC2000 

"BASTZOTO 
3AS02020 

-aASC2030 
8ASC2040 
8ASC2050 
8ASC2C60- 


C  THEN  3313 
33C7  K=K*1 

"T3T9     IF     (X-P)     <=    0    THEN    33CT 
3311     IF     (K-H)>    C 
3313     IF     (I-K)=    0 


THEN    2259 
THEN    2229 


3315    IF     <I-K)>    G     THEN    3359 
2317    FCR    *4=1    TC    M 
3319    TsgJLTtMJ 


BASC2070 
8ASC2080 

3AS0210O 
BASC2U0 

-8A-^C2-t-2«— 
3ASC2130 

_9A_S_e2i4Q_ 
8ASC2150 
BAS02160 

_8A5-C2-W0- 
3ASC2180 
3ASC2190 

BAse-2^ee— 

BAS02210 
_3A_SC22  20_ 
3ASC22J0 
BAS0224O 
3AS02250 
3ASC2260 
8ASC2270 

3AS02290 
_3ASC2300__ 
3ATC2TI0 
8ASC2320 
aA3C23-30— 


T3  21    £( 

23  2  3    8( 
a3Z-5_N4 


I  ,H4)  = 
K,M4)  = 
XT    »4 


6(K,M4) 

T 


1327    13=13+1 
2329    11=1+1 
3321     IF     I1>H 
3333    FCR    M4=I 
3335     IF    e(f4,  _ 
3337    T=B(M4,I)/'e(I,IJ 
2329    FCR    8=1    TC    C 
-3341    8  tfr4,8)  =  E  t>>4,3)-8(  I  ,3  )»T 


THEN    3345 

1     TO    M 

I)=C     THEN    3344 


3343  NE 

3344  NE 
3347  0= 
2349  FC 
3351    C= 


XT  8 
XT  H4 

"t 

R  1  =  1 
C»B(  I, 


TC  M 
JJ 


~JTf3  NEXT  I 
3355  C=( (-1)*EXP( I3)*0 
-3-357  RETURN 


!359  D  = 
3361  RE 
3399  I? 
24C1  Nl 
34C3  FC 
34C5  Dl 
34 C 7  NE 
-34C9  P  = 


C 
TURN 

="t 

=  N+1 
R  1  =  1 
TT=TT(T 
XT  I 
C  =  «  =  Q 


BASC2340 
3ASC2350 

"8ASC2360 
3ASC2370 

3A$_C2  3  8G_ 
3A"SC2390" 
8ASC24G0 

-BA-SC2  440- 
3AS02420 
3ASC2430 
3ASC2440 
BASC2450 
8ASC2460 
8ASC2470 
8ASC2480 

-3A-S42490- 
8AS02500 
3ASC2510 

3ASC2530 
_8A_SC2540_ 
BASC2^50 
8ASC2560 
-BA-SC2  540- 
BASC2580 
8ASC2590 

BASC2610 
BASC2620 

~8AS02630 
3ASC2640 

-8A-S-C26£0- 


TC    M 

r 


3411  IF 
0  THEN 
3413  M 
3415  V ( 
3417  FC 

~j4mrr 

3421  NE 
-3-423-&C 
3425  IF 
2427  IF 
3429  R- 


C(  1) 
3425 
-M-l 
N1)=U( 
R     1  =  1 

tj=t7<t 

XT   i 
TC    341 


NO  I    EwUAL     10    ZERO    THEN    3425 


M)  =  C 
TC    M 


♦XT 


(Nl-1 
(Nl-2 

'Lll)/C12> 


)=C    THEN    3615 

)    NCT    ECUAl    TC    ZERO    THEN    3433 


2421  GC 
2433  IF 
"  3~4~3  5~T= 
2437    w= 

-3439    GC 


TC  354 
IN  1-3 
BT277C 
CUJ/C 
TC    356 


J    NCT    ECUAt    TC    ZERO    THEN    3441 

m 

(3) 


3441    IF     (A8S(C(N1-10/D(N1)-A8S(C{2)/0< 1) )     >=    0    THEN    3437 
34<3     I3=-I3 
3445    f  =  M/2 


344  7  FC 
3449    N2 

~145~I— F= 
3453    C( 

-J-453    CI 


R     1=1 
=  M  +  1- 
CfN2)  " 
N2)1C( 
K2-UCX 


TC    M 

I 


52 


3455  C  (  I  )  =  P 
3457  NEXT  I 
345*9  If  8  NOT 


3461 

JAAA. 

2465 
3467 
14 


P  =  C 
GC  TC 


■£€  IAL  TO  2ESC  THEN  34*5~ 


3465 


3465  IF 


F~=PTC~ 
C  =  l/C 

fl=C  7t-£N  3473 


3471  R=T/R 
3473  E=5E-1C 

3475  JfNl)=CtM) 

3477  W(M)*C(M) 

3 479  J  ( IS  1  +  1)  =  MM»1  )=Q 

3465  FCR  L=l  TC  1000 
3467  FCR  11=1  TC  Ni 


3469  I-NT-U 
3451  J  <  1  J  =  C(  I)*P»J(I+1) 
3453  W(I  )*J(I  )»R*h(  1  +  1) 
N£>T 


3AS02660 
BASC2670 
SASC^6-frQ- 
3ASC2650 
BASC27C0 

BASC2720 
-&ASOZZ3a 
BASC2740 
BAS02750 
8A50276O- 
8ASC2770 
BASC2780 
8Ar02T<30 
3ASC28C0 
8ASC2810 
BAS0282Q 


3455 
3457 

3459 


N~E>T  II 

IF  (AB5U(  l)/0tl  )- 

IF  h.2)  NOT  eCUAL 


E><  =  0  THEN  3549 
IC  ZERO  THEN  3505 


8ASC2830 

8ASC2840 

"3"AS"C28yO" 

8ASC2860 

3AS028JQ- 


35C1  R  =  R*1 
35C3  GCTC  35C7 
?5C5  R-B-J11)/H<2) 


8AS02880 
BASC2850 

3AS02^t3- 


35C7 
35C5 


FCR  11=1 
I=M-I  1 


TC  N3 


3511  JJI  )»C(TT=T!*JU*1)-C*J(I*2J 
3513  l»lI)  =  J(II-f«WI  +  ll-e»MI+2J 
3515    NEXT    I  1 


3ASC2910 

BAS02920 

3ASC2930" 

BASC2540 

8ASC29S0 


3517    IF    C12)    NOT    ECUAL 

3515    IF     (A8St J(2)/0(l) ) 

3521     IF     (AGS(J(2)/Q (1)  * 


TO    ZERO    ThEN    3523 
_E)     >    THEN     3527 
E)     <»    Tt-EN    3525 


BASC2960 
BASC2970 

3ASC2980 


3523  IF 
3525  IF 
r5T7 


(AeSI  J(2)/DU)  )■ 
(A3_U(  1)/0(1)  )- 


3525 

-3-5. 


L_I 


rr="¥(2)-j(2J 

0  =  w(3)*E>F(2)-Cl*Vi 
£ 0  NOT  ECUAL  TC 


E)  >  0  THEN  3527 
E  )  <*  0  THEN  3567 


(4) 

ZSfiC  THEN  3539 


BASC2990 
BASC30QO 
3ASC301~0~ 

3ASC3020 
-3ASa303CL- 


3533 

3535 


P  =  P-2. 

C  =  C*  tC-fl ) 


3537  GCTC  3543 

3539  P=P+ <J(2)*W (3) -J(  1  )*W(4))/C 

3541  Q=d»(-J.2)»Cl+J (  l)«h(2)  )/0 


BASC3040 

3ASC3050 

3AS0306O- 

BASC3070 

BASC3C80 

BA5 03090 

BASC310O 

BASQ3-110 


"3543" 
3545 


7qYXT~L" 
E=E*10 


147  GCTC  3465 


3549  N1  =  M-1 
3551  V(M)=0 


3»-  53 

155 


ir  n  >-  c 

U(N1)»1/R 
TC  3561 


Tl CN  3559 


BAS03120 
8ASC3130 
-3A-5  0314O 


557  GCTC 


3559  U(M)=R" 

3561  FCR  1=1  TC  M 

T.f  3  DJ  I  )  =  .UI+U . 


3ASC3150 

BASC3160 
"3"AiCJl/CT 

BASC3180 
-8ASQ3190- 


3564 
3565 


NEXT 
CTC 


I 

3425 


3567  Nl-M-2 

3565  IF  13  >=  0  TFEN  3577 

3571  U 1=1  /C 

3573  P1  =  P7TC*T. 

3575  GCTC  3581 

-J_5_7_7__CO=C 

3575  Pl=P/2 


8AS032CO 
BASC3210 

-*A-SQ3-220— 
8ASC3230 
3ASC3240 

~3~ATC3~2"5Cr 
BASC326C 

-aAS0327O-- 
8ASC3280 
3ASC3250 

-3AS£33€G- 
3ASC3310 
3AS03320 

"BASC3330"" 
3ASC3340 

_BA<i-_33_5JL- 
3ASC3360 
3AS03370 

-4A-5&3^80— 
3ASC3350 
3AS03400 

"BA-SCT410— 
3ASC3420 

_ilA.S.C3  43jQ_  . 
3ASC3440 
8AS03450 
3ASC3460 


3561 

— 3^-33- 

3565 

3567 


F=P1*EXP<2)-C1 
IF  F  >=  J  THEN 


3-5^ 


UlM  +  1  )  =  L(M)=-P1 
V  (N1>1)  =  SUR (-F J 


. .  ervi  nt  r=-  vtnt+  rr 

3551  GCTC  36C7 
3_593_  _____  2X-  _N.CJLtC.LAL 


3555 
3557 

36C1 

36C3 

36C5 
_26_ 

3613 
3615 
bGCQ 

5CC1 
5CC2 

5CC6 
5CC8 


-TC  iSEG  ThEN  3595L 


U(M*1)=-SCP  (F) 
GCTC  36C1 
iM+l)  =  -(Pl/ABS(Pl)  )MABS{  P1)*SJR(F)) 


V(M  +  i)  =  V<M)=0 
U(M)  =  C1/U(NH-1) 

FCR  1*1  TC-PO 

C(  I  >  =  J(I*2) 
MEXT  I 


GCTC  3425 
RETURN 

*6TLRN 

PRINT 
PRINT 


5010 
5012 


-Ftir-T=rTC- 

PRINT 
NEXT  I 


ThE  ROOTS  ARE 


REAL 


INAGINARY" 


l:uil),vm  ) 


PEINT 

RETLRN 


3ASC3470 
BASC3480 
"SASC3490— 
3ASC35C0 
8AS03510 
JASC3520- 
3ASC3530 
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